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Abstract 

The 3+1 spacetime formulation of general relativity is used to investigate the transverse waves prop- 
agating in a plasma influenced by the gravitational field of Reissner-Nordstrom black hole, as explained 
in an earlier paper, to take account of relativistic effects due to the event horizon. Here, a local approx- 
imation is used to investigate the one-dimensional radial propagation of longitudinal waves. We derive 
the dispersion relation for these waves and solve it numerically for the wave number k. 
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1 Introduction 



In recent years there has been renewed interest in investigating plasmas in the black hole environment. 
The Coulomb potential of charge particles due to coupling is much stronger than the gravitational po- 
tential and are often neglected in the Newtonian approximation. But the mean gravitational held for 
certain astronomical objects like galactic nuclei or black holes may be strong and the observation of 
magnetic fields indicates that a combination of general relativity and plasma physics at least on the level 
of a fluid description is appropriate. A successful study of the waves and emissions from plasmas falling 
into a black hole will be of great value in aiding the observational identification of black hole candidates. 
For this reason, plasma physics in the vicinity of a black hole has become a subject of great interest in 
astrophysics. In the immediate neighborhood of a black hole general relativity applies. It is therefore of 
interest to formulate plasma physics problems in the context of general relativity. 

In the preceding paper pQ, a local approximation has been used to investigate the transverse wave 
propagation in the relativistic two-fluid plasmas surrounding the Reissner-Nordstrdm (RN) black hole 
event horizon. The present paper is concerned with the investigation of longitudinal waves propagation 
in this environment. A general introduction to plasma physics in the presence of the gravitational field of 
black hole has been presented in the preceding paper I, based on the 3+1 formulation of Thorne, Price, and 
MacDonald (TPM) [3] 0] , given detailed in the membrane paradigm [5. , which provides a foundation 
for formulation of a general relativistic set of plasma physics equations in the strong gravitational field 
of black hole. Actually the 3 + 1 approach was originally developed in 1962 by Arnowitt, Deser, and 
Misner [6 to study the quantization of the gravitational held. Since then, their formulation has most 
been applied in studying numerical relativity [7]. 

Exploiting the 3+1 formalism, a lot of researches have been carried out. The study of plasma wave 
in the presence of strong gravitational fields using the 3 + 1 approach is still in its early stages. Zhang 
[8j |9] has considered the care of ideal magneto hydrodynamics waves near a Kerr black hole, accreting 
for the effects of the holes angular momentum but ignoring the effects due to the black hole horizon. 
Holcomb and Tajima [10] , Holcomb 11 , and Dettmann et. al. [12] have considered some properties 
of wave propagation in a Friedmann universe. Daniel and Tajima 13] have studied the physics of high 
frequency electromagnetic waves in a strong Schwarzschild plasma. 

The oscillations of a perfect relativistic plasma, as well as of nonrelativistic plasma, can be divided into 
longitudinal (electrostatic) and electromagnetic (transverse) oscillations. Buzzi, Hines, and Treumann 
(BHT) [TU [15] described a general relativistic version of two-fluid formulation of plasma physics and 
developed a linearized treatment of plasma waves in analogy with the special relativistic formulation of 
Sakai and Kawata (SK) [16 , to investigate the nature of the transverse and longitudinal waves near the 
horizon of the Schwarzschild black hole |14| 115) . In this paper we apply linearized two- fluid equations of 
BHT to investigate longitudinal waves propagating in the plasma close to the event horizon of Reissner- 
Nordtsrdm (RN) black hole. The RN spacetime is the Schwarzschild spacetime generalization with charge 
q and mass M. It is the unique, asymptotically flat, spherically symmetric solution of the Einstein- 
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Maxwell equations that may be analytically extended to an electrovacuum solution representing a black 
hole for < \q\ < M. The extrime RN black hole (i.e., when q = M) is distinguished by its coldness 
(vanishing Hawking temperature) and its supersymmetry. It occupies a special position among the 
solutions to Einstein or Einstein-Maxwell equations because of its complete stability with respect to 
both classical and quantum process permitting its interpretation as a soliton |17[ 118] . The extremal RN 
space is also special in admitting supersymmetry in the contest of N = 2 supergravity [HI [20l ETJ [22l [23] . 
Recently, the quantization of longitudinal electric waves in plasma are studied in |24j and new longitudinal 
waves are found to exist for these types of quantum plasmas [25] . So, our study on longitudinal electric 
waves around a charged black hole is important and meaningful. 

This paper is organized as follows. In section 2, we summarized the nonlinear two-fluid equations ex- 
pressing continuity and conservation of energy and momentum in Reissner-Nordstrdm black hole space- 
time. For zero gravitational field these equations reduce to the corresponding special relativistic ex- 
pressions. In section 3, the two-fluid equations for longitudinal waves have been derived by considering 
one-dimensional wave propagation in the radial z (Rindler coordinate system) direction. We linearize 
these equations for wave propagation in section 4 by giving a small perturbation to fields and fluid param- 
eters. In section 5, we express the derivatives of the unperturbed quantities with respect to z and used 
the local or mean-field approximation to obtain numerical solutions for the wave dispersion relations. We 
describe the dispersion relation for the longitudinal waves in section 6 and give the numerical procedure 
for determining the roots of the dispersion relation. In section 7, we present the numerical solutions for 
the wave number k. In section 8, we present the numerical solutions for the wave number k. Finally, in 
section 9, we present our remarks. We use units G = c = fc^ = I . 



As mentioned in paper [T] , the two-fluid plasma considered here are either electron-positron or electron- 
ion. The equation of continuity and Maxwell's equations for each of the fluid of species s with velocity 
v s , number density n s , and the relativistic Lorentz factor j s in Reissner-Nordstrdm spacetime are given 



2 Two-fluid Equations in RN Spacetime 



by 



— (j s n s ) + V • {a-f s n s v s ) = 0, 



(1) 



and 



V B 



(2) 



V • E 

dB 

~dt 
<9E 

~dt 



(3) 



V x (aE), 



(4) 



V x (qB) - 47raJ, 



(5) 



with the charge q s and the charge and current densities are defined for each fluid species by 





(6) 



S 



S 
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where s is 1 for electrons and 2 for positrons (or ions). The presence of lapse function a defined by 

ot(r) = J = i(r-r+)*(r-r_)*. (7) 

in all the equations signifies the general relativistic effect around a Reissner-Nordstrom black hole. Here, 
r± = M ± y/ M 2 — q 2 . All quantities are measured by a local fiducial observer (FIDO) as discussed in 
paper pQ. The spacetime metric of our coordinates is 



ds 2 = -a 2 dt 2 + -\dr 2 + r 2 {d0 2 + sin 2 <9# 2 ). (8) 

OT 



The gravitational acceleration felt by a FIDO is defined again as 



a = -Vina = -L ( - i-\ e f , (9) 

The energy and momentum conservation equations derived from Maxwell's Eqs. Q and ([5]) coupling 
with each single perfect fluid of species s to the electromagnetic field as given in paper pQ are 

-?rPs - -^Allies + P.)] - V • [ 7s 2 (e s + P s )v s ] + ls q s n s E ■ v s + 2 7s 2 (e s + F s )a ■ v s = 0, (10) 
a at a at 

and 

7 s 2 (e s + P.) + v s • v s + VP S - 7 s g s ^s(E + v s x B) 

+v s (^/ s q s n s E ■ v s + ^-P«J + 7s 2 (e s + P s )[v s (v s • a) - a] =0. (11) 

Here, £ s is the internal energy density and P s is the fluid pressure. If, now, one sets a = 1 so that 
the acceleration goes to zero, these equations reduce to the corresponding special relativistic equations 
as given by SK [16]. Although these equations are valid in a FIDO frame, the transformation from the 
FIDO frame to the comoving (fluid) frame involves a boost velocity, which is simply the freefall velocity 
onto the black hole, given by 

w ff = (l-a 2 )i (12) 

Then the relativistic Lorentz factor 7boost = (1 — W I) -1 ^ 2 = 

The Rindler coordinate system, in which space is locally Cartesian, provides a good approximation 
to the Reissner-Nordstrom metric near the event horizon in the form 

ds 2 = -a 2 dt 2 + dx 2 + dy 2 + dz 2 , (13) 

where 

x = r + (9-~), y = r+tp, z = 2r+a 2 . (14) 

The standard lapse function in Rindler coordinates becomes a = z/2r + , where r + is the event horizon 
of the black hole. 
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3 Longitudinal Two-fluid Equations 

In this section, only a summary of the two-fluid equations relevant for the longitudinal waves will be 
restated. For longitudinal waves it is more convenient to work from a combination of the equation of 
continuity (TTJ), Poisson equation @, and the conservation of momentum equation (|11[) . We can separated 
the longitudinal part of these equations by considering one-dimensional wave propagating in the radial 
z direction. Assuming v SXl v sy and v sz = u s as the velocity components along a;, y and z direction, the 
radial component e z of the equation of continuity ([T]) and the Poisson equation ([3]) can be written as 

d d 

^r(7sn s ) + — (aj s n s u s ) = 0, (15) 
Ot oz 

and 

dE z 

= 47r(gini7i + 92^272)- (16) 

oz 

The two-fluid equations for longitudinal waves can be separated by considering a new complex transverse 
and longitudinal fields and velocities by introducing the following complex variables 

v s (z,t) = v sx (z,t) + iv sy (z,t), v ag (z,t) =u s (z,t), 

B{z,t) = B x (z,t) + iB y (z,t), E(z,t)=E x (z,t) + iE y (z,t). (17) 

These give 

v sx B y ~ v sy B x = 7}(v s B* - v*B), 

v sx E y - v S y E x = l -{v s E* - v*E), (18) 

where the * denotes the complex conjugate. Using the above complex variables, the longitudinal part of 
the momentum conservation equation (lll[) can be written as 

Ps ~Jjr~ = qsTlsls + 5 ( v s B * ~ v *s B )^j + (1 _ u l)Psa 

1 dP \ dP 

g s n s7s E-v s + -±\-—± (19) 

a ot J oz 

where E ■ v s = \ {Ev* + E*v s ) + E z u s and p s is the total energy density defined by 

Ps = j 2 s (e s + P s ) = l 2 s{m s n s + Y g P s ) (20) 

with r g = 7 9 /(7 ff — 1). The gas constant 7 9 take the value 4/3 for T —> 00 and 5/3 for T —> 0. The 
ion temperature profile is closely adiabatic and it approaches 10 12 K near the horizon [26]. Far from the 
(event) horizon electron (positron) temperatures are essentially equal to the ion temperatures, but closer 
to the horizon the electrons are progressively cooled to about 10 8 — 10 9 K by mechanisms like multiple 
Compton scattering and synchrotron radiation. The equation of state can be expressed by using the 
conservation of entropy as 



Where the time derivative in Eqs. dTB]) and (|18|) is defined as 

D_ = (}_d_ i 

Dt \ a dt s dz 



D fid d , 

u s — . (22) 
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4 Linearization 



We linearize the above set of longitudinal two-fluid equations by considering a small perturbation. We 
introduce the quantities 



u s (z,t) = uos(z) + Su s (z,i), v s (z,t) — 5v s (z,t), 

n s (z,t) = n 0s (z) + 6n s (z,t), P s (z,t) = P 0s (z) + SP s (z,t), 

p s (z,t) = p Qs {z)+S Ps (z,t), E(z,t) = 6E(z,t), 

B z (z,t) = B (z) + SB z (z,t), B(z,t) = SB{z,t). 



(23) 



Here, magnetic field has been chosen to lie along the radial ej direction. The relativistic Lorentz factor 
is also linearized such that j s — 7q s + S"/ s , where 



70s 



1 - u os) " > S ls = 7cL u 0s • Su s 



(24) 



Neglecting the product of perturbation terms the conservation of entropy, Eq. (f2Tj) is linearized to 

dP s = bn s . 

Also the total energy density Eq. (|20|) is linearized to 



0p s = H 6n s + 2uos"f sPOsOU s , 

nos V Pos J 

where p Qs = ^ s (m s n 0s + T g P 0s ). 

The continuity Eq. (|15l) . and Poisson's Eq. (fTo]) are linearized to obtain 



(25) 



(26) 



d_ 

di 



+™o s 7 0s 



9 u 0s 
dz 2r_) 

9 9 f 
uos— + a 



d_ | J. 
+ 27+ 



/ f dn 0s 



dt dz 2r + \nos dz 



du, 



5tt s = 0, 



(27) 



and 



dz 



47re(n 27o2 - "•oi7oi) + 47re (702^2 - 701^1) + 47re(n 2Wo27o2^2 - «oiWoi7oV«i)- (28) 



In similar fashion, the longitudinal part of the momentum conservation equation, Eq. (|19[) is linearized 
to give 



du, 



O.s 



+"7o s 



POsJOs 



SE Z 



1 7o 2 s 7 g ^o s ( d d 

7o s ^Os POs V ^ ^ 



73-Pos f 1 ^-Pos 1 dn 0s \ ,f 1+ ihlgPo 



p Qs \P 0s dz n 0s dz 



2 du 0s 



1 



Pos J \ 2r 4 



u 0s ct- 



dz p 0s dz 7o s 2r-+ 



0. (29) 
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5 Local Approximation 

In the paper pQ, the unperturbed radial velocity, magnetic field, number density, and pressure near the 
event horizon for each species s interms of freefall velocity are given by 



uos(z) 

Pos(z) 

so that their derivatives become 



B (z) = B+vUz), 
n 0s (z) = ^n +s Vfi(z), 



(30) 



du Qs _ 


a 1 


dBo _ _ 


4a B 


dz 


2r+ Uff ' 


dz 


2r+ vl ' 


dn 0s _ 


3a n 0s 


dP 0s _ 


3a jgPas 


dz 


2r+ vl ' 


dz 


2r+ vl 



(31) 



where 



c = 



1 



r -(l- r -± 

r + \ r 



and vs = [1 - a 2 (z)} 1/2 . 

As mansion in paper pQ, here we also applying the same restriction on a local scale for which the 
distance from the horizon does not vary significantly and use a local (or mean-field) approximation for 
the lapse function and hence for the equilibrium fields and fluid quantities. Since the plasma is situated 
relatively close to the horizon, i.e., a 2 < 1, we choose a sufficiently small range in z so that a does not vary 
much. In the local approximation for a, a ~ a is valid within a particular layer. Hence, the unperturbed 
fields and fluid quantities, which are functions of a, take on their corresponding "mean- field" values for a 
given ao- 

Then, using the local approximation for a, the derivatives of the equilibrium quantities can be eval- 
uated at each layer for a given a of the form 



DP, 



du 0s 



dz 
dP 



dz 



Dn 0s 



DB n 



dn 0s 



dz 

dB 

dz 



(32) 



The coefficients in Eqs. (|2~T|) . (|2"5|) . and (|2"5|) then become constants within each layer with respect to a 
(and therefore z as well). So it is possible to Fourier transform the equations with respect to z, using 
plane-wave-type solutions for the perturbations of the form ~ e »(fc«-wt) f or each uq layer. 



6 Dispersion Relation 

When Fourier transformed, the continuity Eq. (|2"T|) with the help of Eq. (|3"2")l becomes 



SE 



a k - uju 0s - ia (Dn 0s /n 0s + 3u 0a j$ s Du 0a ) - i/2r + 
a ku 0s - w - i {u 0s /2r + + a ^ s Du Qs ) 



(33) 
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The longitudinal part of momentum equation, Eq. (|29[) when Fourier transformed, gives 



(a ku 0s -ui- ia 7o s ( 1 + u q s ) Du 0s) Su s + 



ia q s n 0s 

A)s70s 



+ 



2 (a fc - wu 0s ) - ia 7 0s — -5— £>"o s £>"o s - 1 H - 

POs POs \"0s nos J \ POs 



a ^o s u OB Du 0s + 



1 

2r~ 



<5n s = 0. 



Finally, the Poisson's equation (|28|) becomes 

ifc(5^ = 4ire(~f 02 6n 2 - 7oi<5«i) + 47re(n 02 wo27o2^ u 2 - noiUoiToi^i)- 



(34) 



(35) 



Using Eqs. (|33p . ([3"4")l , and ([35]), the dispersion relation for the longitudinal wave modes may be put in 
the form 



ao^ ~ 77- i«o7o s f — Dn 0s + 2jn s u 0s Du Qs 

LT± \ riQ s 



ia 7o s (l + u 2 0s )Du 0s \ ( a ku 0s - u) - ia j^ s Du 0s - 



J Ts 



i / 1 2 

a k - ujuqs - iota Dn 0s + 3"/ 0s u 0s Du 0: 

2r+ \n 0s 



1 



1 



a k - uju as - ia Q [ -^-DP 0s - -^-DP 0s ) I 1 + -7^ I I h aoJo s uo s Du 0l 

-njs -fos / v Ts \ z j \ir. 



(36) 



where v\ s = 2j g ~f 2 s Po s / pg s is the thermal velocity, u cs — ejo s no s Bo/ po s the cyclotron frequency, and 
ujp S = 47re 2 7o s ng s /pos the plasma frequency of the fluid's of species s. 



Note that, the plasma frequency is also frame independent, that is, independent of 7o s since the 7q s 
factor in the numerator cancel out the 7q s factor involved the energy density, po s , in the denominator. 
The same is also true for the thermal velocity, for which the 7g s factor cancels. The above dispersion 
relation is valid for either electron-positron or electron-ion plasma in that it makes no assumption as to 
the mass, number density, or temperature of each species. In this respect it is completely general. If one 
considers the equivalent case to that of SK [16] in Eq. (|36]l for an electron-positron plasma, in which the 
two fluids have the same velocity uq , the same equilibrium density rto, and at the same temperature To, 
one obtains 

"'Ml I / 

— a k - ia 7n I 

x < [a ku Q —w — ia 7o(l + ul)Du ] I a ku — to - ia "foDu 



a k - ia o7o ( Z~ Dn v + 27oito^'Mi 



2r+ 



a^k — luuq 



2r, 



ia ( — Dn Q + S^uqDuq 
n 



2i 



O.s 



a k - ujuq - ia ( — Dp o - — DP 0s ) - -5- ( 1 + — ] ( — + a j u Du Q 



2r 4 



(37) 
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In the limit of zero gravity, i.e., a — > 1 and with the velocities uq — > 0, Lorentz factor 70 — > 1, the derivative 
of unperturbed quantities Duq = Diiq = DPq = DBq = 0, it reduced to uj 2 — [2uip + k 2r y g Po/po), which 
recovered the SK [TB] result. The only difference being that 7 g was set to unity in the SK work. 

7 Numerical Solution of Longitudinal Modes 



The dispersion relation given in Eq. (|36[) is complicated enough, even in the simplest cases for the 
electron-positron plasma where both species are assumed to have the same equilibrium parameters, and 
an analytical solution is cumbersome and unprofitable. We solve numerically the dispersion relation in 
order to determine all the physically meaningful modes for the transverse waves. We put the equations 
in the form of a matrix equation as follows: 



(A-kl)X = 0. 



(38) 



The eigenvalue is chosen to be the wave number fc, the eigenvector X is given by the relevant set of 
perturbations, and / is the identity matrix. 

We need to write the perturbation equations in an appropriate form. We introduce the following set 
of dimensionless variables: 

1 



Ul 



k = — , 



r ~ Su s _ Sv s Sn s 

ou s = , v s = , dn s = , 

Uqs u 0s n 0s 

6E Z 



6B ~ SE 
6B = -5-, E=—, 6E Z 
-£>o -do 



B 



(39) 



The derivatives of the equilibrium quantities given in the equation 05] are also expressed in terms of k + 
so that 

«ofc+ (4Q) 



dz 



and 



Du 0s = - 
Dh 0s = - 



3aok + 



DB = 3 — B , 



-n 0s 



DPos 



3a k + ~f g 



(41) 



For an electron-positron plasma, ui p i = lo v i and uj c i — lo C 2, and for electron- ion plasma ui p = ^Jio p \u v i 
uj c = ^/o; c iw C 2, and u>l s = (2uip S + lo 2 s ). Here, w» is defined as 



W» = UJp 



\f{ui v \UJ v -i). 



(42) 



The dimensionless eigenvector for the Longitudinal set of equations is 



X 



Longitudinal 



Sill 

5u 2 
Shi 
6n 2 
SE, 



(43) 
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Using Eqs. (|39[) and (|4Tj) , the Eqs. (|33|). (|34|) . and (|35|) can be written in the following dimensionless 
form 



1 



7_ ~_ ~ 1 "o s w | 1 - -TT I - i — 



767 <, - uff.,/2 2u 0s 7o s 



-<5n. 



i(g s /e)w c 



(44) 



k5n s 



*0s 

2 2 
U 0s70s 



u 0s ^ 1 



<4 s /2 Vm 0s 7o 



i27o s /c H 



• fc +«o / 2 

1 !5— (u, 



Os 



'Os 



/2) * 



i(g s /e)u; c 



w p (ug s - wf.,/2) 



(45) 



2 2 2 2 

• 2 2 W pl r~ -22 W p2 c~ , • W pl r~ • W p2 

m oi7oi — — du i _ m 027o2 — — ' )w 2 + i — - — oni - l- 



(46) 



These are the equations in the required form to be used as input to Eq. (|36|) . Using Eqs. (|44)) - (|461) . the 
Eq. (|38|) then can be written as 

(A - kI)X = 0. (47) 

The numerical solution is carried out using MATLAB as performed in paper [T] with C, 2 — 1.05 and 
q 2 /M 2 = 0.2. The limiting horizon values for the electron-positron plasma and electron-ion plasma are 
taken respectively as 

n +s = 10 18 cm- 3 , T +s = 10 10 K, B+ = 3 x 10 6 G, and lg = -, 

n +1 = 10 18 cm- 3 , T +1 = 10 10 K, n +2 = lO^cm" 3 , T +2 = 10 12 K, 
with 7 5 = |, M = 5Mq , g 2 = 0.2Af 2 . 

8 Results 

8.1 Electron-positron Plasma 

1. Longitudinal Low Frequency Mode 

The longitudinal modes for electron-positron plasma are split into high and low frequency modes with 
domain as Q < \/2 and Co > \/2. For low frequency our study show five modes to exist, while for 
Schwarzschild case four low frequency modes are found |15j . These are in contrast with the special 
relativistic case investigated by SK |16j where only one high frequency mode was found to exist for 
electron-positron plasma. The first low frequency mode shown in Fig. 1 is physical and shows growth for 
«o < ott but is damped for ao > a*, where at denotes the transonic radius occurs at about at ~ 0.82 for 
the case considered here. This means that energy is drained from the wave rather than being fed into it 
by the gravitational field. The second low frequency mode is similar to the mode shown in Fig. 1. The 
third low frequency mode shown in Fig. 2 is also physical for ao < at- This mode is odd in that there 
appears to be some interplay between frequency and distance from the horizon ao, which split this mode 
into two distinct regions. One is purely growth region for which Im(k) >> Re(k) and the other which is 



10 



Figure 1: Longitudinal low frequency mode showing both damping and growth for electron-positron plasma. 




Figure 2: This mode shows some interplay between frequency and distance from the horizon ao for 
electron-positron plasma. 

also a growth region but where Im(k) « Re(k). The fourth low frequency mode is opposite in sense to 
the previous mode and is not shown here. The fifth and final low frequency mode has a little difference 
from the mode shown in Fig. 2 and it is also true for this mode that there appears some interplay 
between frequency and distance from the horizon ao. It is evident that the growth and damping rates 
are independent of the frequency, but depended only on the radial distance from the black hole horizon 
through a - 

2. Longitudinal High Frequency Mode 

In the high frequency domain there exist five modes. The first two modes arc equivalent to the modes 
shown in Fig. 1 and show growth for ao < at and damped for ao > at- Again, energy is drained from 
the wave rather than being fed into it by the gravitational field. The third mode shown in Fig. 3 is a 
new mode which shows damped corresponding to a < a t and growth for a > a t . The other two modes 
(one is depicted in Fig. 4) also show growth for a < a t and damped above a t . 
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Figure 3: Longitudinal high frequency growing and damped mode for electron-positron plasma. 




Figure 4: Longitudinal high frequency damping and growing mode for electron-positron plasma. 

8.2 Electron-ion Plasma 

1. Longitudinal Low Frequency Mode 

For electron-ion plasma the longitudinal modes are split into high and low frequency domain as u> < 1 
for low frequency, and u> > 1 for high frequency. As for electron-positron case, here we have found five 
low frequency modes to exist. Since the temperature of electron and ion near the horizon are different so 
their transonic radius is same for-electron positron case occurs at a t \ ~ 0.82 . The transonic radius for 
ion occurs at a t 2 >^ 0.99 and can not properly be shown for low frequency modes. The first two modes 
are similar to the modes shown in Fig. 4 for electron-positron plasma and both show growth for a < a t \ 
and damped for a > a t \. This is interesting again that energy is being fed into the wave between the 
transonic radii but is drained from the wave closed to the horizon. The third mode illustrated in Fig. 5 
is damped for < cto < 1. The remaining two modes are equivalent to one another and show growth 
for a < a t i and for a t \ < a < a t 2 and only one of these modes is shown in Fig. 6. It is clear that 
the growth and damping rates are frequency independent as for electron-positron plasma, but depended 
only on the distance from the black hole horizon through ciq. 

2. Longitudinal High Frequency Mode 
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100, 




Figure 5: Longitudinal low frequency damped mode for electron-ion plasma. 




Figure 6: Longitudinal low frequency growth mode for electron-ion plasma. 

Four high frequency modes are found to exist as shown in Figs. 7 and 8. The first mode not shown 
here is similar to the first high frequency mode shown in Fig. 1 for electron-positron plasma. The second 
mode shown in Fig. 7 is also growing for ao < ctti arid damped for ao > an- There appear some interplay 
between frequency and lapse function ao- Thus near the transonic radius, it appears that energy is no 
longer fed into wave mode by the gravitational field but begins to be drained from the waves. Fig. 8 
shows two modes, one of which shown in bottom left is growing for a < a t \ and damped for a > a n , 
while the bottom right one is damped for a < a t 2 and growth for a > a t2 . The transonic radius for ions 
is clearly shown here. The instability (damped and growing) of longitudinal wave modes has occurred at 
each transonic radius and is interesting. 

9 Concluding Remarks 

The prime concern of this study has been exclusively the investigation, within the local approximation, 
of longitudinal waves in two-fluid plasma situated around the event horizon of the Reissner-Nordstrom 
black hole. The transonic radius begins to play a significant role for the longitudinal waves. From Eqs. 
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Figure 7: Longitudinal high frequency mode both damping and growth for electron-ion plasma. 



(|44]) and ([45|) it is clear that a singularity take place for each fluid at the point for which the freefall 
velocity equals the half of the thermal velocity, Mq s = \v^ s . The position of transonic radius of each fluid 
mainly dependent on their limiting horizon temperature and determines the temperature at any given 
radius. Near the transonic radius damping and growth occur for all the low and high frequency modes. 

Using a local approximation the dispersion relation for longitudinal waves has been derived. In 
the limit of zero gravity this result reduced to the special relativistic result obtain by SK [16] for the 
longitudinal waves. One interesting point concerning the longitudinal waves in the electron-positron 
plasma is that, unlike the result found by SK |16j for which only one high frequency mode exist. Here 
there is no such restriction on the frequency and five low frequency modes are found for both the electron- 
positron and electron-ion plasmas. Different modes become physical (Re(k) > 0) at the boundary define 
by the transonic radius of each fluid. It is true for the majority of the modes except for two intriguing 
low frequency modes for electron-positron plasma that some complicated interplay between frequency 
and redial distance appears to determine the regions in the u) — olq plane for which the mode is physical. 

The presence of magnetic monopole charge in the RN hole and the characteristic of the extremal 
RN hole draw attention of the physicists. In view of these reasons, our study of longitudinal waves 
propagation in a relativistic two-fluid plasma in the environment close to the event horizon of the RN 
black hole is interesting. The result we obtained reduces to that of the Schwarzschild black hole [15] when 
q = 0. Our result can be specialized for the extreme RN hole by choosing M 2 = q 2 . In view of these 
reasons, our study presented in this paper is thus well motivated. 
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